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Abstract. We study a class of nonlinear diffusion equations whose model is the classical 
porous media equation on domains Q C R*, N > 3, with homogeneous Neumann boundary 
conditions. Firstly we improve some known results in such model case, both as concerns sharp 
L q °-L°° regularizing properties of the evolution for short time and as concerns sharp long time 
CN) ' asymptotics in the sense of L°° convergence of solutions to their mean value. The generality 

of the discussion allows to consider, almost at the same time, also weighted versions of the 
above equation provided an appropriate weighted Sobolev inequality is required to hold. In 
fact, we show that the validity of a slightly weaker functional inequality is equivalent to the 
<**] . validity of a suitable L qo -L°° bound for solutions to the associated weighted porous media 

^ ' equation. The long time asymptotic analysis relies as well on the assumed weighted Sobolev 

inequality only, and allows to prove uniform convergence to the mean value, with the rate 
£Nj ■ predicted by linearization, in such generality. This fact was not known even for the explicit 

0\| classes of weights previously considered in the literature. 
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1. Introduction 



< 

& • The aim of this note is to provide some new L qo -L°° regularity and asymptotic estimates for 

solutions to nonlinear diffusion equations whose model is the following Porous Media Equation 
with homogeneous Neumann Boundary conditions: 



u t = A (u m ) in n x (0, oo) 

= on 8Q x (0,oo) , (1.1) 
^(•,0) = tio(-) in 



where ft is a bounded regular (say C 1 ) domain of M. N , m > 1, iV > 3 and, as usual, for y ^ we 
define y' m : = |y| m_1 y when dealing with non-necessarily positive solutions. Even in this case, in 
. which we shall prove results which will turn out to be sharp both for short and long time, such 

estimates will improve those obtained in the pioneering paper [T] and, for larger classes of data, 
in [5] (see O [2] for the Neumann problem in the case in which £1 has infinite measure) . Actually, 
almost all of these results will also work for the Weighted Porous Media Equation introduced 
in |13j (see |17[ \TE\ \TT\ \12\ [26l [271 E] an d references quoted therein for the one- weight case) : 



p u u t = div (p^ V(u m )) in Q x (0, oo) 

Pm^? = ° on^x(0,oo). (1.2) 

k u(-,0) = ito(-) in fl 

Indeed, the only relevant assumption for our method to work is that the weights p u ,p^ are 
strictly positive in f2 (but may be degenerate or singular at the boundary), sufficiently regular 
and such that the Sobolev-type inequality (assuming throughout the paper that z/(f2) < oo) 

\\v-v\\ 2a . v <C s \\Vv\\ 2 .^ V^W 1 ' 2 ^;^) (1.3) 

holds true for suitable Cs > 0, a > 1, where v is the mean value of v with respect to the 
measure dz^ = p u dx (we also set d/x = /f^dx) and || • || p: ^ will denote, for any pG [1, oo), the L p 
norm w.r.t. a given measure A on f2. In fact the validity of the weaker inequality 

\M*r,u < C 's (ll V«|| 2yi + IMIiJ Vv e W^(n; p) (1.4) 
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will be sufficient for most of our purposes, and actually will turn out to be equivalent to the 
following bounds for the L°° norm of solutions that we shall prove: 

ll«(*)lloo < K Ir^-v^-v Unoll^^- 1 ' + \\u \\ qo J Vi > 0, 

K being a constant which depends only on go € [1, oo), m, Cs, o~ and 

We shall discuss these issues in greater detail in Section [5] One should anyway compare these 
kind of results with the ones of |13] , which involved problems similar to (|1.2p under the weaker 
requirement that the Poincare inequality 

\\V -V\\ 2 ;v < C P\\^ V hv W 6 Wl ' 2 {&\ V,fi) (1.5) 

holds true, so that min [S (— Aw") \ {0}] > 1/Cj> > 0, where S (— Ajy) denotes the L 2 (£l;v) 
spectrum of (minus) the weighted Laplacian — Aj^(v) = — p~ l div (p M Vi>) with homogeneous 
Neumann boundary conditions. While of course the results of [13] continue to hold in the present 
situation, stronger ones will be shown to be valid and, in particular, solutions corresponding to 
L qa data become instantaneously bounded with quantitative bounds on the L°° norm, a fact 
which need not be true under the sole condition (|1.5p . as shown in [13] . 

As for the long time asymptotics, we discuss first the non- weighted case and show uniform 
convergence of solutions to the mean value of the initial datum with a sharp rate. Again 
this bounds improve on the results of [5]. In the weighted case, we shall prove that uniform 
convergence to the (weighted) mean value of the initial datum occurs, for general L 1 data, as 
a consequence of the sole validity of the Sobolev-type inequality (|1.3p : the proof of this fact is 
completely different from the one known in the non-weighted context and requires a much more 
delicate functional argument. In particular, if the weighted mean is not zero, the rate we give 
is exactly the one predicted by linearization. The existing results (see |18| HU [10]) show, in the 
explicit class of weights considered there, only pointwise or local uniform convergence; however 
we stress that the Sobolev-type inequality we require need not always hold in their setting. 
As already mentioned, we comment that uniform convergence does not necessarily occur if one 
assumes only the validity of (jl.5p . as shown in |13] by explicit counterexamples. 

For a thorough study of smoothing and decay properties of solutions to large classes of 
nonlinear evolution equations on see the monograph |29| . whilst for other work specifically 
concerned with the connection between functional inequalities and asymptotic properties of 
solutions to weighted porous media equations we refer the reader to (EMU EI], remarking however 
that L go -L°° regularization properties of the evolutions considered are not addressed there. The 
pioneering papers |18| . |14] show local uniform convergence of solutions to their mean in some 
explicit one dimensional, single weight case but do not deal with rates of convergence. For the 
Neumann problem in the case in which f2 has infinite measure (where convergence to zero is 
considered) see e.g. [3] [2] . 

The proofs of existence and uniqueness of weak solutions to problem (jl.ip can be found in 
[1], |30| . A proof of existence of solutions (and of uniqueness in the class of energy solutions) to 
(|1.2p has been given in |13| . under certain regularity conditions on the weights: in particular the 
problem is well-posed provided the weights are sufficiently smooth and locally bounded away 
from zero in (while they can be singular or degenerate at the boundary), a fact that we shall 
assume hereafter. 

Our first main result, Theorem 13.21 gives a bound on the L°° norm of solutions to (jl.ip for 
short times and L qo (go > 1) initial data. This bound improves considerably Theorem 1.1 of 
[5] (see the discussion in Remark 13. 3p and it is sharp in the sense that it captures the exact 
explosion rates of the well-known Barenblatt solutions [30] for short time. Bounds for large 
times are given in Theorem 14.11 for zero-mean data, they being sharp (as already noticed by 
[1] for L°° data) and improving the corresponding results of [5]: in particular, an absolute 
bound for the L°° norm follows. The case of data having non-zero mean is also considered in 
Theorem 14.31 where we show that convergence to the mean value occurs with an exponential 
rate matching the one predicted by linearization. Once again we improve the results of [5] and 
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in Proposition 14.51 we show the sharpness of our bounds by producing a class of data for which 
a matching lower bound for the rate of convergence holds. 

The corresponding results for solutions to (|1.2p are given in Section [5J see in particular The- 
orems [5^21 E21 ETO] and (5TTTJ Besides proving such results, which are parallel to those obtained 
for the non- weighted case, we show in Theorem 15.31 and in Corollary 15.41 that the validity of 
a short time L 90 -L°° regularizing effect for solutions to (jl.2p implies in turn the validity of a 
suitable Sobolev-type inequality like (jl.4p . thus giving a converse to the corresponding regu- 
larization theorem and hence showing its optimality, in some sense. Uniform convergence to 
the (non- zero) weighted mean of the initial datum is shown in Theorem 15.101 whereas Theorem 
15.111 gives an explicit rate of convergence. 

Finally, Section 15.11 provides a concise list of explicit classes of weights for which suitable 
Sobolev inequalities hold. We remark that almost none of the corresponding nonlinear diffusions 
seems to have been studied in the existing literature. 

2. The concept of solution 

For the reader's convenience we begin recalling from [13] (see also [T] and [30] for the non- 
weighted case) the concept of solution we are going to consider. We write it in the weighted 
case only, but clearly the definition applies for the non- weighted case as well just by setting 
Pu, Pfi = 1- Recall that du = p u dx and d/i = p M dx. 

Definition 2.1. A function 

u G L 2 ((0, T); L 2 ($7; u)) : V(u m ) G L 2 ((0, T); [L 2 (Q; p)] N ) VT > 

is a weak solution of (|1.2p with initial datum uo G L 2 (f2; v) if it satisfies 

n«(x,%(x,t)di/dt = - [ tt (x)r?(x,0)dz^ + f [ V(u m )(x, t) ■ Vr/(x, t) dfidt 
a Jn Jo Jn 

Vt? G W^ 2 ((0, T); L 2 (fl; v)) : V? ? G L 2 ((0, T); [L 2 (ft; fj)] N ) , r](T) = . 

Uniqueness of solutions does not hold in general in the weighted context, even for bounded 
data and solutions, however it does indeed under the additional condition that 

u e L m+1 {(0,T);L m+1 (n;v)) , 

see |13[ Prop. 3.10]. We call such a solution, if any, a weak energy solution of the equation 
considered, according to a common terminology. In fact Theorem 3.12 of |13| gives a proof 
of existence of the weak energy solution for initial data uq G L rn+l (ft; v) , under some further 
regularity assumptions on the weights, for example Pu,Pfi > and 

For more general data uq G L l (VL;v) the correct extension of the concept of solution is given 
in \12>\ Sect. 3] (the discussion being similar to the one first introduced in [301 Sect. 6.1] for 
the non- weighted problem). Such solutions are called limit solutions, since they are naturally 
obtained as limits of energy solutions, that is by approximating no with a sequence of data 
Uo n G L mJrl (Q;u) and exploiting the fundamental L 1 -comparison principle, namely the fact 
that if u and v are the solutions to (|1.2p with initial data respectively uq and vq then 

\\{u(t)-v(t)) + \\ lrv < ||(« -«o)+|| 1; * V *>°- (2- 1 ) 

Thanks to (|2. 1 1) . in particular, the sequence of energy solutions {u n } corresponding to {uQ n } 
is Cauchy in L°°((0, oo); L 1 (fi; v)). For more details about the mentioned existence results of 
energy solutions, we refer the reader to |13| Th. 3.12] or to the original |30| Th. 11.2] in the 
non-weighted case. 

Throughout the whole paper, when referring to "the solution" to the equation at hand, we 
shall mean without further comment the unique weak energy or limit solution constructed as 
above. Also, we shall often make use of two fundamental properties of such solutions. In first 
place, the mean value of the initial datum is preserved along the evolution: see |13| Prop. 3.13] 
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for a proof. In second place, the inequality ||tt(i)|| p;l / < ||w(s)||p ;l , holds true for allp € [1, oo] and 
all t > s, as a consequence of the energy estimates of |13^ Sect. 3]; we shall refer to the latter 
property as "non-expansivity" of the norms. For the corresponding results in the non-weighted 
case see [3D] Th. 11.10]. 

In the following sections, when dealing with (|l.ip . Q will always be a bounded C 1 domain of 
M. N . As for the weighted case (|1.2|) regularity of $7 in principle is not needed, although it may 
be hidden in the assumed validity of the corresponding weighted Sobolev inequalities. 

3. Short time L qo -L°° bounds 
We start recalling a useful numerical lemma given in )13] , 
Lemma 3.1. Given a, f3 & (0, 1), with a > (3, there exists a constant c = c(a, (3) > such that 
x -a y l-a + X -P y l-P + y < c ( Qj p)(x-<*y l - a + y) Vx, I/6I+. 



Proof. We need to show that 



R(x,y) 



x a y 1 a + y 



is bounded in ]R + x ]R + by a constant which depends only on a and j3, and in order to do that 
one finds explicitly the zeros of R x (-,y) for any given y. □ 

Now we prove one of the main results of this paper which, as we shall remark later, is a sharp 
improvement of the L qo -L°° regularity estimate first provided by Theorem 1.1 of [5]. 

Theorem 3.2. Let u be the solution of (|l.ip corresponding to an initial datum uq € L qa (£l), 
with qo £ [1, oo). The following estimate holds: 



IK^IL < k (0,0+^-1) |K||; 90+JV(m - 1) + KIM vt > o, (3.i) 

where K is a constant which depends on m, Cs, N , and can be taken to be independent of 

go- 

Proof. We shall proceed by means of a classical Moser iterative technique. Firstly we consider 
an initial datum no € L°°(f2): the fact that the estimate we shall obtain will not depend on 
1 1 ^0 1 |oo w iU allow us to extend it to general L q ° data thanks to a well-known argument that we 
shall recall at the end of this proof. 

Given t > 0, let us consider the sequence of time steps % = t (1 — 2~ k ). Clearly, to = and 
too = t. Also, let {pk} be an increasing sequence of positive numbers such that po = qo and 
Poo = oo, which we shall explicitly define later. For the moment, we assume in addition that 
qo G (1, oo) fl [m — 1, oo) (afterwards we shall be able to remove this hypothesis). Multiplying 
(jl.ip by ii Pfc_1 and integrating in f2 x (tfc,tfc+i) (and neglecting ||w(tfc+i)||pfe) gives 



4(fffc ~ l)Pkm 
{Pk+m- l) 2 



k+l 



p k +m-l 

V « 2 x,s 



dxds<\\u(t k )\\ p p k k . (3.2) 



With no loss of generality, suppose |fi| = 1. In order to suitably handle the left hand side of 
(|3.2p . it is convenient to notice that the validity of the Sobolev inequality 

\\v-v\\ 2N < CaWVvL Vuew 1,2 (n) 

JV-2 

implies the validity of the inequality 

5^NlLr-7^NI?<l|V«||^ v B e^(fl), (3.3) 



POROUS MEDIA EQUATIONS WITH HOMOGENEOUS NEUMANN BOUNDARY CONDITIONS 5 
where a = N/(N - 2). Upon applying fl3J3]) to the function u (Pk+™-i)/2 in (jjj^ we get: 



2(p k - l)p k m /"* 



)Pkm /" <fe + 1 r 

Since go > m ~~ 1 an d Pfc i s increasing, we can control ||it||( Pt ,+m_i)/2 with ||tt|| Pfc . By that and 
by the non-expansivity of the norms, we deduce: 

fa - - l) W m 2 _ fc y < || M(<t )| C + 7 ^ 2-- « HtlJIia^- ■ 

UsKPk + m — 1) ^sKPk + 771 — 1) 

(3.4) 

provided pfc+i = <r(pfc + m — 1). Now wg assume ||iao||oo — This hypothesis, again together 
with the non-expansivity of the norms, ensures that 

so that (|3.4p reads 

Pk+i a 2 (r>, -A- <m — 1 "l 2 

n^ + i)ii P ; +1 < t k -D P J 2 r + 2 ii«c**)iis • ( 3 - 5 ) 

Clearly there exists a suitable constant D = D(qo, m, Cg, N) such that (|3.5|) simplifies to 



C|(p fe + m - 1) 

<K**)II5 + t? 



4(p fc - l)p fc m Z"*^ 1 | Li/ „M, Pfc +m-i 



lh(4 + i)lU +1 ^d^^ + i)^ iiu(t fc )n p ; fc+i . 

Setting Uf. = \\u(tk)\\ Pk , it is straightforward to check that the sequence {U^} satisfies 

a k + 2 -(k+2)a + k + l <j k + 2 -a g k + 1 

u k+1 < d pfe+i^- 1 )" (r 1 + 1) ffc+i^- 1 ' u qo Pfc +! . (3.6) 

Also, one can verify that p k = (qo — A)a k + A, with A = — — Ti) = y(l — Ti). Letting 
fc — )■ oo in (|3.6p . we get (from now on D = D(qo, m, C$, N) will denote a generic constant which 
may differ from line to line): 

N ?TO 

\\u(t)\\oo = lim Ht)\L < hminf U k+1 < D{t~ l + 1) ^o+^-D ||„ ||- . (3>7) 

fc->00 + 1 fc— >0O ^0 

Actually, (j3.7[) is not an L qo -L°° regularity estimate. Indeed, recall that it has been obtained for 
initial data uq such that ||uo||oo = 1- Through a simple time scaling argument we can deduce 
from it an estimate for general L°° data. That is, given a solution u(-,t) to (jl.ip corresponding 
to the initial datum uq, it is straightforward to check that u(-,t) = ju(-, A 1 "" 1 t) is the solution 
of the same equation corresponding to the initial datum no = j u o- Choosing A = ||uo||oo an d 
applying (|3.7p to u(t) we conclude that 

N 2q 

h(t)\L = hoimhoWZ' 1 1) < D(t~ l + ||uoC" 1 ) a * +JV(ra - 1) KIIJ D+JV(m " 1) • (3-8) 

Estimate f|3 . 8[> . as such, is not of great interest. However, it is possible to reduce its dependence 
on 1 1 Mo | |oo i n the following way. First, let us rewrite it as 

that is, setting — 



2q +N(m-l) ' 

IK*)lloo < D (r^r||n ||; o - e + WML • (3.9) 

Now it is convenient to exploit the semigroup property, by shifting the time origin from to 
t/2, and the non-expansivity of the L qo norm. This leads us to 

< d (r^r + \\ u (t/2)C WuoW]; 9 ) ; 
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applying to ||u(t/2)|| oo we obtain 



IK*)IL < d [t'^\\u^ 6 + t-^\\u Q \\\f + h4t Hlf 2 

it is then clear that proceeding in this way along n steps one arrives at 

K^lloo < D(n, •) (r^ \\uoW 1 - 9 + ... + r^r WuoW]; 9 " + \\u \C hoW]; 9 ") ■ (3.10) 

In order to remove the dependence on the L°° norm, we need a suitable L 90 -L°° regularity 
estimate. To this end, one can reason in the following way. Suppose ||tio||<j = 1- From (|3.4|) , 
setting again U k = ||^(ifc)||p fe , we have: 

u k+1 PJ ^ < ° 2 j Pk + 7 " 1)2 2 fc r 1 ^ + 2CV* 4 "- 1 • (3.H) 
(pfc -l)p k m 

Consider a solution to (|3.1ip . namely a sequence {14} such that Vo = C/o = 1 and 

V k+1 ^ = ^ + 7" 1)2 2* r + 2V k ^ ; 

it is easy to verify (by induction) that U k < V k and V k > 1. Proceeding similarly to the case 
H M o||oo = 1 (now the leading term on the right hand side is V k Pk+m ~ 1 ) we deduce that for a 
suitable constant D = D(qo,m,Cs, N) the sequence {V k } also satisfies the following recursive 
inequality: 

jc + l a 

V k+ i < D Pk + 1 {t' 1 + l)"w-i V k . (3.12) 

Solving (|3.12p and exploiting the fact that p k > o~ k one gets (recall that D may vary from 
line to line) 

(fc + 2)(l- CT ) CT - fc - 1 - t T- fc - 1 + CT a-a-k 
V k + l<D 90("-D a (t~ l + 1)«0("- D , 

and passing to the limit as k — > oo, 

= lim ||u(i)L < liminf f/ fc < liminf V fc < L> (t _1 + 1)^ . 

k— ¥oo * k k— ¥oo k— >oo 

Again, by means of a time scaling argument (now with A = ||ito||go)> we nave that the final 
L qo -L°° regularity estimate is 

iV(m-l) 



jV 1 



which reads, upon setting 5 = N(m — l)/2qo, 

Ht)\\co < D (*~^ WMl; 5 + IKIlJ • (3.13) 

Of course the exponents involved in f|3. 13[) are not very satisfactory, in particular when go is 
small with respect to N. However, combining this estimate together with (|3.10p one can obtain 
a much better L^-L 00 regularity result. Indeed, by means of the usual shift of the time origin 
to t/2 in (|3.10p and thanks to (|3.13p evaluated at time t/2, we have: 

IK*)l| oo <^(-,n)(r^Klt + --- + r^ I IKII^ +t—\\u \\ 1 - 5e +\Wo\\ qo ) ; 

now we pick n great enough so that 50 n < 9, and apply iteratively Lemma [3. II with x = i 1 /^™ -1 ), 
V = \\ u o\\q , at = and (3 = 59 n at the first step and then f3 = 6 J along j = n . . .2. We thus get 

ggg 



< D l t - 2qo+N(m -i) ii^n^o+^-i) + || Uo || \ (3 w) 

which is valid, as stated at the beginning of the proof, for go £ (1) °°) H [m — 1, oo). 

We now extend the above estimates to the case of a general go > 1- In fact, for notational 
simplicity, we shall give the complete proof only for go = 1 . We shall exploit a technique similar 
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to the one used in the proof of Corollary 8.1 from [24]. That is, first of all consider the analogue 
of estimate (j3. 14|) in the time interval [t/2,t], for qo = m: 



< D 



f \ 2m+iV(m-l) 



\u(t/2)\\ 



2m+N(m-l} 



+ ||«(t/2)|| 



using the inequality 
obtain 



< 



< D ( rn-l 



\oo 1 ^ m || • ||i^ m and the non-expansivity of the L 1 norm we 

(m-l)(l-7) 1 — 1 m-1 1 



where 



|«(V2)|loo m KHi™ +KV2)|| C 
N(m - 1) 



Pol 



i 



(3.15) 
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2m + N(m-1) 

In order to handle (j3. 15j) . we can argue as in the previous part of this proof. That is, consider 
first an initial datum uq such that ||uo||oo < 1- This, in particular, implies that ||it(t/2)|| (X3 < 1 
and ||tfco || 1 < 1) so that (|3.15p becomes 



(m-l)(l-7) 1-7 

<D (t m-i+1) ||«(V2)IL m Unollx" 1 , 



(3.16) 



up to absorbing 2" / ^ m ^ into D. It is apparent that iterating estimate ()3. 16j) along k steps one 
gets 



Halloo <2 



1 yph=k-l h ( (m-l)(l- 7 ) * 
m— 1 £sh=0 \ m 



D t 



^ +1 



h=fc-l / (m-l)(l-7) ' 



Eh=K 
h=0 



X 



it (t/2 



'(m-l)(l-7)' 



1-7 y^h = fc-l ( (™,-l)(l-7) ' 



(3.17) 



Polli 



En=K 
h = 



Passing to the limit in (|3.17f> as k — > oo we end up with 

1-7 

\u 



1-7 



< D [ r— 1 i+("-ih ||u ||i 1+Cm " 1)7 + IKo|li 1+(m_1)7 



(3.18) 



for another suitable constant D. Estimate (|3.18p is analogous to (|3,7p . By reasoning likewise 
we can obtain again (|3.10p . now with go = 1 and 9 - 



i+(m-l)7 ' 

< D{n, •)(* "-Mpolli + t m - 1 ||«o||i + Folloo 



iii 
l^olli 



(3.19) 



In order to remove the dependence of the right hand side of (|3.19p on ||tio||ooi it is convenient 
to suppose first that ||ito||i > 1 and look for a suitable L l -L°° regularity estimate. Indeed, 
considering again (|3.15p . it is easy to check that the worst possible case occurs when 1 1 C s ) II oo > 1 
for all s < t: assuming that, clearly ||u(t/2) ||^ 7 ^ m < ||it(i/2)||^ 1 '' m ^ and since also 



n o|li 



< ||tto||] 7 ' m we have that ()3.15p reads 



m — 1 
| m 
I OO 



\Ht)\\oo< D [t m -' + 1 ) IK*/2)IU" 

From (|3.20p . through analogous computations as above we easily get 

/ hit \ 
<D\t rn-l \\uq\\ 1 + ||«0|ll ) , 



(3.20) 



(3.21) 



which holds provided ||ito||i > 1- By means of a time scaling argument (recall how we obtained 
(13~T3D ) dMH) becomes 



* — 1 Kill 



1 — TO7 



-ollij , (3.22) 

which is valid for all uq E L 1 (f]). By choosing n in (|3.19p great enough so that m^9 n < 9 and 
then combining (|3.19p with (|3.22p through the usual t/2-shift argument, we obtain: 

1-7 \ 

(3.23) 



<D[t \\uo\\{ 



l + (m-l)7 



+ INIIi , 
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which is exactly (|3.14p with qo replaced by 1. As previously stated, a similar strategy also works 
when qo € (l,m — 1], provided this latter interval is not empty. Following the above proof it is 
easy to realize that the constant D in the inequality corresponding to (|3.23p when || • || 9o instead 
of || • ||i is considered, is bounded as a function of qo £ [1, m — 1]. Moreover, one can notice 
that the constant D in (|3.14p is locally bounded as a function of qo € (1, oo) n [m — 1, oo) and 
it remains bounded as qo oo. Hence (|3.ip holds for all qo > 1, with a multiplicative constant 
K independent of qo- 

Finally, we are left with removing the hypotheses = 1 and uq € L°°(f2). For the first 
one it is enough to proceed through a spatial scaling argument exactly as explained in the end 
of the proof of Theorem 5.4 of |13j . Namely, if u(x,t) is a solution of (jl.ip in a domain Q of 
measure with initial datum tto(x) then 

u(x,t) = \n\~ N (.™-Vu(\n\^x.,t\ (3.24) 
is also a solution of (jl.ip in the domain = Jl/|Jl|^f of measure 1 with initial datum 

U (x) = |0|~iV(m-l) Uo f|0|^x) . 

Noticed that, one applies (|3.ip to u and goes back to the original solution u by means of (|3.24p 
and the relations 

\\u\\ p = \n\~w^-p \\ u \\ p , c s (o) = |n|3-^-wc s (n) . 

In consequence of that, it all amounts to admit the dependence of the multiplicative constant 
appearing in (|3.ip on too. 

The extension to general L qo data is handled by means of a standard argument. That is, 
given uo 6 L qo (Q), consider a sequence {«on} C L°°(f2) which converges to no in L qo (£l) and 
the corresponding solutions to (jl.ip {u n } and u. From the L 1 -comparison principle we deduce 
that, for any given t > 0, u n (t) — > u{t) in L 1 (il); moreover from (|3.ip we also learn that, up 
to subsequences, {u n (t)} converges in the weak* topology of L°°(£l) to an element w G L°°(Q). 
The identification between u(t) and w is straightforward, and estimate (|3.ip is preserved to the 
limit thanks to the weak* lower semicontinuity of the L°° norm. □ 

Remark 3.3. We stress that estimate (|3.ip in fact improves, with respect to the dependence 
on the time variable and on the parameters qo, m, N, the one first provided by Theorem 1.1 of 
[5] (obtained by means of a Gross differential method), which reads 



\u(t)\\„ < Cr~ \\uo\\l~ a e £oll " ol| iv (m -i)* ) (3.25) 



where C,Eo are suitable constants that depend only on <7o> Cg, N, \Q\ and 

jv- 



1 



<]() 



qo + m — 1 



(3.26) 



Indeed, it is plain that 



N( m — I) 

a> - ^Mf u Vm>l, ViV>2, V 90 >1. 
2^o + N(m - 1) 

Also, note that (|3.25p blows up as t — > oo. Moreover, it only holds for initial data belonging to 
L qo (Q), with qo > IV (m — 1). Instead, estimate (|3.ip does not blow up and holds for all initial 
data which lie in L qo (Q) for any go > 1. 

Remark 3.4. Consider the well-known Barenblatt solutions (recall that we work for > 3) 
. , A („ ,|x| 2 \^ , N A A(m-l) 



t 2 ~< J ' A(m-l) + 2' ' N ' 2mA 
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which solve the Neumann problem (jl.ip for sufficiently short times, provided x = belongs 
to (else one just translates them). By a scaling invariant estimate for solutions to (jl.ip we 
mean a bound of the type 

IK^IL < s(t, Kill) vt > o , Vtto g L\n) 

with S(-, •) such that 

S(t, \\uqWi) = IWohSiWuoW^-H,!) Vt>0, VuoGL 1 ^). 

So S(-, •) is completely determined by the one- variable function S(y, 1). Estimate (|3.ip is clearly 
scaling invariant. Moreover, it is a simple calculation to check that the left and right hand sides 
of (|3.ip . when evaluated on the Barenblatt solutions and considered e.g. in the time interval 
[t, 2t], have the same rate of divergence as t J, 0, and that the L 1 norm of u#(-,i) is preserved 
along the evolution. This means that f)3. 1|) (for go = 1) is sharp in the sense that there cannot 
hold another scaling invariant estimate with a better rate for S(y, 1) as y J, 0. Also, we shall 
see in section [5] that, should the bound 

K^lloc^^^^lKII^ + IKIlJ Vt>0, Vn GL 90 (O) (3.27) 

be valid for some 6 < 2q ^+N(m-i) anc ^ ^ or an ^° — ^ ^ n f ac t> Qo £ [ m i m + 1) would suffice), then 
an embedding of W 1,2 {Q) into L 9 (S7) with q > 2N/(N — 2) would be true as a consequence, a 
fact which is clearly false since N >3. 

Remark 3.5. We point out that we have stated our estimates for all t > rather than for 
almost every t > 0. This remains correct since in |30| Ths. 11.2, 11.3] continuity of solutions 
to (jl.ip in C([0, T]; L 1 (0)) (for all T > 0) was proved, while in |13] (see Section 3) the same 
was shown to hold for solutions to (|1.2p at least when < oo (which here will be the case), 
of course upon replacing L 1 (Q) with L 1 (0;i/). This comment also applies to the results of the 
forthcoming sections. 



4. L°° BOUNDS FOR LARGE TIME 

In this section we shall prove some sharp improvements to the asymptotic estimates provided 
by Corollaries 1.3 and 1.4 of [5]- We remark that both the L^-L 00 regularity estimates (|3.ip 
and (|3.25p give no significant information about the behaviour of the solution u(-,t) as t — > oo. 
So in order to obtain such informations one has to proceed through different techniques, whose 
starting point is the distinction between zero mean and non-zero mean solutions. 

4.1. The case u = 0. The following result, as we shall discuss in greater detail later, provides 
an improvement of the asymptotic estimates given by (5J Cor. 1.3] for zero-mean solutions. 

Theorem 4.1. Let u$ € L qo {Q), with go € [1, oo) anduQ = 0. For the solution u of (jl.ip with 
initial datum uq the following estimate holds: 

N 1 

IK*)lloo<Qi* 2qn vt>o, (4.i) 

f r. . , II ||l-m\ (m-l)(2 90 +JV(m-l)) 
Q2t + Hliollgo J 

where Q\ and Q2 are constants depending on qo, m, SI. In particular, the absolute bound 

IK^IL^Qa* - ^ Vt>0 (4.2) 

holds true. 

Proof. We start proving the assertion for q$> 1. In this case, the bound 

Ht)\\ qo < " -JT (4-3) 

Bt + \\u \\l; m 
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(B = B(qo,m,Q)) was already established in the proof of Corollary 1.3 of [5]. Actually, this 
was done only for qo > 2, however one sees easily that the same result can be obtained also for 

go e (1,2). 

Now, let us consider estimate (|3.ip with the time origin shifted to t/2: 

\\u(t)\\oo < K ^y^^ \\u(t/2)\\^^ + \\u(t/2)\\ q ^ . 

Exploiting (j4.3|) evaluated at time t/2 we remain with 

I _ „ \ 

IKt)IL<* 



t\ 2, +JV(m-l) 1 1 

— + 



9 / / \ ?22_ 

> / TDt , II Ill-m\ (m-l)(2 90 +J* 

V (^2 + IKIU J ^2^H"0|| 30 ; J 

from which (|4.ip follows immediately. The absolute bound (|4.2p is a straightforward conse- 
quence of (|4.ip . Notice that performing the above calculations in full detail one gets that 
Qi — > oo and Q2 — > as qo X 1) thus making it impossible to perform a limit to handle the case 
qo = 1. So in order to prove (|4.ip for qo = 1 we need a different argument. First we notice that 
the validity of such inequality for qo = 2 implies that 

. « 1 1 ^ tq ,, ,,i_<v n N(m — 1) 

ll«(0lloo < Qi* "-MKIIs 70 



To 



4 + iV(m- 1) ' 



hence we have also 

TO 



(f\ m-l 1-7Q t-TQ 

2) ll^/2)|loo 2 ll^oll! 2 • (4.4) 

Using the strategy outlined in the proof of Theorem I3.2[ see in particular the calculations 
following (|3.16p . one can show that the validity of (|4.4p implies in turn the validity of 

IN^II^Qsi-^lluoll^ , uh>- N ( m ~V 



2 + JV(m-l) ' 
Therefore it follows that 

ll«(t)lla < \Wo\\! < QI t-^^WuoW 1 -^ . (4.5) 



Exploiting (j4. 1[) on the time interval [t/2,t] with the choice qo = 2 and ()4.5|) at time t/2 yields 

1 

(Q 2 t 1 + \\u{t/2)\\ 



7Q 

f\ m-l 1 

K*)lloc<Ql ( o) T^< 

\l—rn\ m-l 



< 



t \ ^1 



9 / \ 1-7(1 

Q22+Q3 \2) W U 0\\l 



t\ m-l 1 



9 , \ !~TQ 

TO - <"0 / 1 — m »o (-1 „.)(-] u 0\\ m—1 

(|)^ (02 I + <3 3 "^ (&* Kill )( t) 



15 



1-7| 

Q 2 (|) 2 +Q 3 2 Kill 



Hence (|4.ip follows also for qo = 1 upon redefining the multiplicative constants and using the 
fact that for any given 6 G (0, 1) one has 2 8 ' 1 (a 5 + b 5 ) < (a + bf < a 8 + b s for all a, 6 > 
(<5 = 1 — Uq/2 to our purposes). □ 
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Remark 4.2. In [5j Cor. 1.3] the following estimate was proved: 

IK*)lloo< v*>i, (4.6) 

(^-i) + |KIIJr) m_1 

for suitable constants B, C > 0, where a € (0, 1) is defined as in (|3.26p . Clearly (|4.6p is weaker 
than (|4.ip . since it provides us with a slower decay rate for ||it(t)||oo as t oo. Actually, a 
rate of order t~ 1 /( m ~ 1 ) for t large was already known to hold for bounded initial data since the 
pioneering work [1] (see Theorems 3.1 and 4.1 there), even though the estimate given there still 
depended in a nontrivial way on the L°° norm of the initial datum and in particular was not 
extendible to data not belonging to L°°{Q). The sharpness of such rate follows by considering 
the separable variable solutions constructed in [TJ Sect. 2]. 
Finally, notice that (|4.ip for small t is equivalent to (|3.ip . 

In the statement of Theorem 14.11 we pointed out that the constants Qi, Q2 depend on the 
whole 0, (in which we then absorb their dependence on TV and Cs as well). The reason traces 
back to the constant appearing in Lemma 3.2 of [1] (exploited in the proof of [5] Cor. 1.3]), 
from which B in (|4,3p depends, whose existence is established by means of an argument by 
contradiction and so whose relation with significant features of the domain (such as its Poincare 
constant) is completely unknown. We recall that Lemma 3.2 of [1] states that there exists a 
constant Cq > depending on r > 1 and SI such that for all integrable functions £ such that 
£ = and £ r 6 W 1,2 (Q) the following inequality holds: 

un 2 <cbiivo 2 - 

In fact, the result is also true for r > 1/2. In |13[ Lem. 5.8] a generalization of the above lemma 
is given, although still obtained through an argument by contradiction. 



4.2. The case u^O. We now show an improvement of the asymptotic estimate provided by 
[5] Cor. 1.4] for solutions with non-zero mean. When making use of the Poincare inequality 
(|1.5p we shall implicitly assume that the constant Cp is the best one for which such inequality 
holds: in other words 1/Cp = X\, where Ai is the first non-zero eigenvalue of (minus) the 
Laplacian operator with Neumann boundary conditions. 

Theorem 4.3. Let u be the solution of (jl.ip corresponding to an initial datum uq € L 1 (0) 
with «o = u 7^ 0. There exists a constant G = G(m,Cs,N, \Sl\,u) > such that the following 
estimate holds: 



\u(t) 



^lloo<Ge 



-\u[ 



Vt > 1. 



(4.7) 



Proof. We start considering data in L 2 (Q). From the uniform convergence to the mean value 
proved in [5] Cor. 1.4], we deduce in particular that there exists a time tq = tq(u) such that 



inf |it(x,t)| > — Vt > Tq . 
xgn 2 



(4.8) 



Using (|4.8j) in the differential inequality solved by \\u(s) — u\\g (for a fixed g > 1), one gets 



d_ 



us 



-m 



11 



m—l 



V \u(s) 



. e 

U\ 2 



Vs > Tq . 



(4.9) 



In the case g = 2, (|4.9p immediately leads to 



\u{t) 



< e p 



(t-ro) 



\ u i T o) 



Vt > T , 



(4.10) 



upon applying (|1.5p to the right hand side. Now we observe that, setting tf. = TQ+(t— tq)(1 — 2 



and p fcH 



N 
N-2 



Pk (let pq = 2), integrating inequality (|4,9p between tf~ and tk+i with g = Pk 
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and proceeding with similar (but easier) computations as in the proof of Theorem 13 . 2 1 one arrives 
at the following estimate: 



I — i -i m—l 



N 

-1 \ T 



(t-T ) +1 ||u(t )-«|| 2 Vt>T , (4.11) 



for a suitable constant D = D(m,Cs, N, |f2|). By replacing To with ^y 2 - in (|4,lip . exploiting 
(pTTUj) at time ^ (instead of t) and using the non-expansivity of \\u(s) — u\\2 (consequence of 
(|4.9p) we obtain 

IK*) -u\\oo < Ge~ Mt \\u -u\\ 2 Vt>r + 1, (4.12) 
for another suitable positive constants G = G(m,Cs, N, \£l\,u) and M = M(m,Cp, \u\). Let 
us take T\ > To + 1 such that 

^\u\<\u\-Ge~ Mt \\u -u\\ 2 Vt>n. (4.13) 

Since 

|ii(x, s)\ > |u| — |it(x, s) — u\ > \u\ — Ge~ Ms \\uo — u\\ 2 Vs > t\ , for a.e. x E , 

we have: 
d 



ds 



\u(s) -uf 2 = -2m [ |it(x, s)| m_1 |Vu(x,s)| 2 dx < 



<-2m(|n|-Ge- Ms ||n -n|| 2 ) m 1 ||V«(s)||^ < (4.14) 

2?n , —Ms n n l n / \ n2 lj ^ 

< — — J- (p| — G e ||tto — u\\ 2 ) \\u{s) — u\\ 2 Vs > t\ . 

By integrating (|4.14p from ri to i (and again exploiting the non-expansivity of \\u(s) — u\\ 2 ) we 
get: 



\\u(t)-u\\^ 2 <\\u -u\\ 2 e c p tiK ' . (4.15) 

Setting e = Ge~ Ms \\uq — u\\ 2 , thanks to the elementary inequalities (recall (|4. 13j) ) 

(\u\ - e) m ~ l > \u\ m ~ l - (m-l)\u\ rn - 2 e Vm > 2 , 

(|«| -e) m ~ X > Inl" 1 " 1 - (m- 1)(|«| - e) m " 2 e > Inj™" 1 - 2 2 ~ m (m - l)\u\ m - 2 e Vm G (1,2), 
and to a straightforward computation of the time integral, estimate (|4,15p reads 

||«(t) -«|| 2 < Ce Vt>n, (4.16) 

being G another suitable constant, different from the one above but depending on the same 
quantities (we can absorb the Poincare constant Cp into C5). Combining (|4.12p with the time 
origin shifted to t — tq — 1 and (|4.16|) evaluated at time t — tq — 1 (in place of t) and setting 
r = tq + T\ + 1 we finally arrive at 



-^lup- 1 * 

loo 



IK*) -till,*, < Ge c p Vt>r, (4.17) 

again for some G = G(m, Cs, N, [fi|, u). The passage to general L 1 data and the fact that one 
can take, up to a different constant G, r = 1 (or any other fixed r > 0) in (|4.17p follow from 
the L^-L 00 regularizing effect. □ 

Remark 4.4. Estimate (|4.7p improves the one given in Cor. 1.4] with respect to time 
dependence. Indeed in [5] the authors provided a decay rate of order 

- il-")™ \u\ m -H 

GC-± om-l 1 1 
e c p 2 

for \\u(t) — Uo||oo (a is still defined by (|3.26|) ). which is clearly slower than the one obtained 
here. However, we remark that an exponential rate of the type 

~e^r\u\ m - l t 

G £ e °> Ve G (0, 1) 
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could have also been obtained by exploiting the uniform convergence to the mean value proved 
in [5] together with [TJ Th. 3.3], although the constant G £ would diverge as e t 1. The new 
contribution of our Theorem 14.31 lies in the fact that one can actually choose e = 1. Notice that 
writing u(t) = u + f(t) and linearizing around the constant solution u one formally gets 

the equation ft = m\u\ m Af, so that the rate given in Theorem 14.31 coincides with the one 
predicted by linearization. 

In the following we shall prove that estimate (|4.7|) is sharp by providing data for which a 
lower bound matching the upper bound in f)4.7j) holds. The result we shall give depends on the 
well-known fact that, under our assumptions, the operator — Aj\f (non- weighted) possesses a 
sequence of eigenvalues {A n } (possibly repeated) and a corresponding sequence of eigenfunctions 
{4> n } which form an orthonormal basis for L 2 (J7). 

Proposition 4.5. Let £1 be a bounded C°° domain. Consider the first non-zero eigenvalue X\ of 
(minus) the Laplacian with homogeneous Neumann boundary conditions, and the corresponding 
set of eigenfunctions (possibly consisting of one single element) normalized in L 2 (Q). Choose 
any of such eigenfunctions and denote it by ifj\ . Given uq ^ 0, let u be the solution of (|1.1|) 
corresponding to the initial datum uq = uo + ciipi, with \c\\ sufficiently small. Then there exists 
a constant H = H(m, $7, c\) > such that the following estimate holds: 

||u(*) -Silo,, > He c p Vt>0. (4.18) 

Proof. Let us consider first the solution to (jl.ip corresponding to a regular initial datum no 
such that uq = 1 and ||tto — l||oo < $ < 1 ; with 5 to be chosen sufficiently small afterwards. 
Since \\u(t) — l\\<x> is non-increasing (indeed \\u(t) — l\\ g is non-increasing for all g > 1, and so 
for g = oo as well) we have that |u(x,t)| > 1 — 5 for all x € and t > 0. Hence by inserting 
such bound into the first line of (|4,14p and solving similarly the resulting differential inequality 
we obtain the following estimate: 

||/(i)|| 2 <ll/o|| 2 e- Aim(1 - 5)m ~^ Vt>0, (4.19) 

where we set /(x, t) = u(x,t) — 1 and /o(x) = uq(x) — 1. From standard quasilinear theory 
(see |22j) wc have that / € C°°(f2 x (0, oo)). Also, straightforward computations show that / 
solves the following differential equation: 

f t = mAf + F, (4.20) 

up to defining F as 

F = m(m - 1)(1 + f) m ~ 2 \Vf\ 2 + m [(1 + f) m ~ l - l] A/ . 
If we choose as initial datum /q (x) = c\ tp\ (x) , under the condition 



|ci| < -n-r-n— , (4.21) 



n oo 



we see that (|4.19p becomes 



2 < |ci| e - Xim ^ m * Vt > . (4.22) 

Define as usual |/|^o(q) := ||/||oo and, for any multi-index r\ = (rji, . . . ,rj^f), the quantity 
\r]\ = Tlx + . . . + t]n and the seminorms 

l/lc*(iD : =™^IWIIoo. k . 

From the uniform parabolicity of the equation at hand we infer the existence of a constant 
Q = Q(5) > such that 

\f(t)\ ck(n) <Q Vt>0,V£:eN, (4.23) 

where Q can be taken to be independent of c\ subject to (|4.2ip . We now recall the generalized 
interpolation inequalities 

N+jp p(k-j) 

\9\ cm < C jAp \g\%» \\g\\r kp , (4.24) 
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valid for all integers k > j > and real p > 1 (see \23\ p. 130] or, for a short review, [6l App. 
3]). Combining (|4.22p . (|4.23p and (|4.24|) we deduce that the following bounds hold: 



ll-e e -Aim(l-5) m - 1 (l-e)t 



for/ = 0,1,2, Vt>0, VeG (0,1), 



being Q'(5, e) > another constant independent of c\ subject to (|4.2ip . From the very definition 
of F and from the fact that the L°° norm of f(t, •) is non-increasing we have that, for suitable 
constants Qq = Qo(S), Q" = Q"(5,e), 



\F 



< 



Qo(\f(t)\ 2 c i m + \\f(t)\L \f(t)\c<R) 



< 



< Q"|ci| 2(1 " e) e -2A im (i-5)-Mi-e) t \/i > , Ve £ (0, 1) . 



(4.25) 



We want to study the asymptotic behaviour of a±(t) = (f(t), iPi)l 2 (Q)i that is the first Fourier 
coefficient of / along the evolution. In order to do that, let us multiply equation (|4.20p by 
and integrate on f2, so as to obtain the following differential equation for a\{t): 

di(t) = -\\ma-]_(t) + (F(t),ipi) L 2(ty . 

Duhamel principle entails that a\ must satisfy the following nonlinear integral equation (recall 
that by construction «i(0) = c\): 

rt 

Ai ms 



ai(t) 



Ximt 



Cl + 



(F(s),ipi) L 2 {n) ds 



(4.26) 



If we choose 5 and e sufficiently small, namely such that 2(1 — S) m 1 (1 — e) > 1, thanks to 
(|4.25p we easily see that the time integral in (|4,26p can be bounded in the following way: 



Aims 



<B\ci 



,2(l-e) 



where B is a suitable constant independent of c\ and t > 0. Pick now |ci| small enough so that 
| ci | > B \ci\ 2 ^~ e \ this being possible since 2(1 — e) > 1. Under such bound on \c\\ we then 



~ e 



deduce from (|4.26p that the exact decay rate of ct\(t) is e Aim *, in the sense that |«i(i)| 
as t — > oo. Hence ()4.18p holds true since there exists a suitable constant K > such that 



Ximt 



W(t) - 1|L = \\f(t)\\oo > ll/(i)ll 2 > l^l" 5 l«i(*)l > Ke 



-Ximt 



Vt > 0. 



The case of initial data with non-zero mean value uq 7^ 1 can be brought back to the case uq = 1 
by means of a standard time scaling argument. □ 



5. The weighted case 

As already mentioned, the main results provided in the previous sections can be extended to 
the case of the weighted porous media equation with Neumann boundary conditions (jl.2p , dealt 
with in detail in |13| . under the sole hypotheses that the weights p u ,Pfj, satisfy Sobolev-type 
inequalities like (j 1 . 3 1) or the weaker (|l,4p . Since the proofs of the corresponding theorems basi- 
cally would be the same as for the non-weighted case, we shall only give statements with short 
comments, except in the cases of Theorem 15.31 Theorem 15.101 and related lemmas, which re- 
quire new arguments. Before doing that, we briefly recall some notation and essential definitions 
concerning weighted Sobolev spaces. 

Let f2 C M N be a domain and let v and \i be two measures defined on it, both absolutely 
continuous with respect to the Lebesgue measure. We indicate as p u and the corresponding 
weights (or densities), which will always be assumed to be strictly positive almost everywhere. 
For all p S [1, 00) we introduce the Banach space L p (f2; v) of equivalence classes of Lebesgue 
measurable functions / such that 

11/11^ = \\ft v dv=l |/PV,dx<oc. 
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The same applies for the measure p. According to |20] . we define the weighted Sobolev space 
W 1,p (Cl;u, p) as the space of all (equivalence classes of) functions v E W,' (0) such that 

\\ V \\ P p;u,„ = W V W P p;u+ II Vv fp;n < °° ■ 

Without further assumptions on p v and p^ in general W ,p (f2; i/, p) would not be complete. 

Definition 5.1. For all p E (1, oo) we denote as B p (£l) the class of all measurable functions g 
such that g > a.e. and 

f^eL} oc (n). 

One can prove [HI Th. 2.1] that if p € (1, oo) and p M E B P (Q) then W^ p {^;v,p) is indeed 
complete. If p = 1 the same result holds true providing that the condition p^ E S p ($7) is 
replaced by p^ 1 € L£.(fi). 

The fact that for any ip E C£°(Tl) the quantity ||y||p;i/u is finite is equivalent (see \20\ Lem. 
4.4]) to the local finiteness of v and p, that is 

pu,Pn,e l} oc {Vl) . (5.1) 

Assuming (|5.ip , one can define the space Wq' p (Q; v, p) as the closure of C™{0) in W la> (Q; v, p). 



Now we are ready to state the weighted counterpart of the L qo -L°° regularity Theorem 13.21 
Hereafter, the assumptions on the weights (inner regularity and boundedness away from zero) 
and the concepts of solution will always be the ones given in Section [2j 

Theorem 5.2. Let z/(f2) < oo and let inequality (jl.4p hold true for some a > 1. Then for 
the solution u of (jl.2p corresponding to an initial datum uq E L qo (£l;v) with qo E [1, oo) the 
following estimate holds: 

( <T (<7-l)<iO \ 

K*)IL < K Ir^^ ll^o|| ( l 7; )90+CT(m ~ 1) + ||n || go; J Vt > 0, (5.2) 

where K is a constant which depends only on m, Cs, o~ and z/(fi). 

The proof of Theorem 15.21 is exactly the same as the one of Theorem 13.21 indeed, going back 
to the latter, one easily sees that in order to repeat it the only relevant assumption is the fact 
that the weaker inequality (|1.4p holds (with a = N/(N — 2) in that case). 

Now we show that the L qo -L°° regularity estimate (|5.2p is in fact the correct one, namely its 
validity implies in turn the validity of a Sobolev embedding for W 1,2 (Q; v, p). 

Theorem 5.3. Let < oo and suppose that there exist a constant K > and two given 

numbers a > 1 and qo E [m,m + 1) such that, for all uq E L qo (Q;i/), the solution u of (|1.2p 
corresponding to the initial datum uq satisfies the following estimate: 

\\u{t)\\oo < K U (CT - 1)90+CT(m - 1) \Wo\\^ )qo+a{m - 1) + \\u4 qo -A Vt > 0. (5.3) 
Then there exists a constant C > such that the functional inequality 

IML;„ < C ' (ll^H^ + IMIfaJ V« E W^ 2 ^l;u,p) (5.4) 

holds as well. 

Proof. We proceed along the lines of the proofs of Theorems 4.3 and 5.6 of [13] , First of 
all, consider a non-negative initial datum uq E L°°(Cl) n W l,2 (VL; u, p) and its corresponding 
solution to (|1.2p u. Exploiting the classical interpolation inequality between the exponents 
qo E [to, to + 1), to + 1, oo and the non-expansivity of the L qo (£l\u) norm, one easily obtains 
that 

m + l-q q 

;iloo m+1 ll u o||^ 



Ht)\\ m+ i:u < IW*)IL m+1 hlS Vt>0; (5.5) 
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combining (|5.5p with estimate fj5.3[) one gets, up to a multiplicative constant that for the sake 
of simplicity from now on we shall always denote as K, 

qQ (2am — in — 1) 



<t)\CX\.<K 



er(m + l — qn ) 

t - ( CT -i), + CT (^-i) i| Un i|( CT - i ^o+ CT (" 1 - 1 ' + ik n ir +i 

II "U ||on:f ~ II u II go;'-' 



Mt > 0. 



(5.6) 



Now we need to use the following fundamental inequality: 



Ht)\ZXl» ~ WMZXiu > -(m + l)t IIVWH^ Vt > 0, (5.7) 



which was proved in |13| (see again Theorems 4.3 and 5.6). From 
that 



and (|5.7p we deduce 



I nm+l ^ 
l M 0|| m +i ;i y — 



< K 



CT(m+l-q ) 
t (ct-I^o+cO™- 1 ) 1 1 1*0 



go(2frm — m — 1) 
(CT-l)9o+°"(m-l) 



+ (m + l)t||V(^)||| M +|| Uo 



im+l 
I go;" 



V* > 0. 



(5. 



Minimizing the right hand side of (j5.8[) with respect to t > one arrives at 



l nm+l ^ ts 



1^0 



<2q (2irm- m— 1) 
I 2am — qQ 

I go;" 



|V«* 



2p-(m+l — <jq) 
I 2crm-(jQ 



+ "0 



im+l 
1 90)1' 



(5.9) 



Proceeding exactly as explained in the proof of Theorem 5.6 of [13] . that is by approximating 



l/ra 




with a sequence of regular functions, one can show that (|5.9p is equivalent to the following 



inequality: 
llf IL... < K 



|Vu|| 2 + \\v\\<m 

m. ' 



Vv > G L°°(n)nW 1,2 (S2;i/,/i), i = - + - — -, 

r q s 



where it is understood that 

m + 1 



(5.10) 



2cr, S 



2am(m + 1 — go) 



m m (m + l)(2<rm — go) 

The fundamental result [H Th. 3.1] ensures that (up to a different multiplicative constant) one 
can put i? = 1 and thus r = 2a in (|5.10p , this leading to the claimed inequality (|5.4p at least for 
regular positive functions; the fact that such inequality holds in the whole space W l,2 (Vt; u, fj,) 
then follows by means of a density argument and by writing v = v+ — V-. □ 

A straightforward consequence of the above Theorems 15.21 and 15.31 is the following, which is 
one of the main results of the paper. 

Corollary 5.4. Let u(Q) < oo. The validity of inequality (jl.4p for some a > 1 is equivalent 
to the validity of the family of estimates (|5.2p for all go € [1, oo). 

Remark 5.5. We shall provide a specific example in which the bound (|5.2[) (for go = 1) is 
sharp in the sense that no better scaling invariant estimate can hold, as in the non-weighted 
case discussed in Remark 13.41 

The setting is one-dimensional, in particular we choose Q = (0, 1). We consider the weights 
Pu{x) = 1, Pfi{x) = ar . The Sobolev inequality (|1.3p associated to such weights is known to hold 
with a = -g^- if /3 £ (1,2) (see Section [5.11 for appropriate references). An explicit calculation 
shows that the Barenblatt-type functions 



UB,p(x,t) 



t -C [ C - k—r 



i 

m— 1 



(5.11) 



where C > is a free (mass) paramenter, £ 



m+l-j8> U * ^ n m(2-/3)(m+l- / 9) ' 

are solutions to the corresponding weighted equation for any t > 0. We consider them for t 
small enough to ensure that the support of Ub a('i~t) is bounded away from the point x = 1. 
It is straightforward to prove that the mass of UB,p(-,t) is preserved in time, so that also the 
Neumann boundary condition at x = is satisfied. One then checks directly that such solutions 



2-/3 and k 



m—l 
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belong to the appropriate functional spaces, so that they are energy solutions (for t > 0) in the 
sense of Definition 12.11 

We notice that (minus) the power of time appearing in the bound fj5.3[) coincides, for q$ = 1 
and a = ^zy, with m+ ^_^ , namely with the power of time £ appearing in (|5.1ip . The fact 
that there cannot hold another scaling invariant estimate with a better rate then follows as in 
Remark 13.41 

We want to stress that Corollary 15.41 has a natural counterpart for the Dirichlet problem 

'ut = p^ 1 div( / 9 M V(n m )) innx(0,oo) 

< u = on dn x (0,oo) , (5.12) 

it(-, 0) = uo(-) in O 

whose well-posedness was also analysed in [13] and whose correct L 9o -L°° regularity estimates 
basically had already been proved in [5j Th. 1.5] (there is no difficulty in extending the proof 
provided therein to our context). Here is the statement: 

Corollary 5.6. Let i/(f2) < oo. The validity of the Sobolev inequality 

M\2*-,u<C S \\Vv\\ 2ill WG< 2 (fi;^) (5.13) 

for some a > 1 and some constant Cs > is equivalent to the validity of the family of estimates 

IKOIL < Kt-Wo+'^-v K||£- 1)?0+CT(m - 1) Vt>0, Vgo€[l,oo), (5.14) 

where u is the solution of (j5.12|) corresponding to a generic initial datum uq G L qo {Vt]v) and 
K > is a suitable constant independent of Uq. 

The finiteness of the measure, in this case, is necessary only in order to prove the converse 
(i.e. the fact that flgjg)) implies (pTTHj) . see p3J Th. 4.3]). 

As for the asymptotic behaviour of solutions, we begin with an analogue of Theorem 14.11 

Theorem 5.7. Let u(p,) < oo and let the Sobolev-type inequality (|1.3p hold true for some a > 1. 
Then if uq G L qo (£l; v), with go G [1, oo), uq = and it is i/ie solution of (|1.2p corresponding to 
the initial datum uq, the following estimate holds: 

imoiioo < < 9i*" (<r - i), ° ;<r(m - i) 1 - vt>o, 

/'n + I II \\l-m\ (m-l)[(cr-l)q +v(m-l)] 

[Q 2 t + \\u \\ qo . u J 
where Q\ and Qi are constants depending on qo, m, £1, v, [i. 

Again, the proof of this result is similar to the one of Theorem 14.11 Actually, in order to be 
able to exploit Lemma 3.2 of [I], a priori one should require that the embedding W ' (fi; f, fJ,) 
L 2 (Q;v) is compact. In principle this is not a major restriction since it is quite common that, 
at least on bounded domains, the validity of a Sobolev-type inequality comes together with 
the compactness of such embedding (see the examples of Section 6 of |13j). Nonetheless, as 
previously mentioned, in |13[ Lem. 5.8] a generalized version of that lemma is provided, which 
does not require compactness. 

As concerns solutions to (jl.2p with non-zero mean, in order to obtain the weighted version 
of the asymptotic estimate given in Theorem 14.31 first we have to prove uniform convergence of 
such solutions to their mean value. Since the proof of Cor. 1.4] (which takes advantage of 
the uniform spatial Holder continuity of solutions to (jl.ip ) in the generality of present context 
does not work, we shall provide a different (functional) argument, which basically only uses the 
assumed validity of the Sobolev-type inequality (|1.3p . To this end, we need two preliminary 
lemmas. 

Lemma 5.8. Given r > 1/2, m > 1 and a fixed constant R > 1, let 

$r, m ( x )= I \y\ r - l \y + i\^ &y Vxe[-R,R]. 
Jo 
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Then there exist two positive constants C\ = C\{m) and C2 = C2(m,R) such that 

1+1 v[^] <\*r,m(x)\<y\x\ r Vxe[-R,R]. (5.15) 

Proof. One has to study the ratio |$ r , im (x)|/|x| 7 *. The bound from above is easily obtainable, 
since 



771 — 1 



|$r,m(s)| <(R+1) — 



yr'dy 







7 71—1 

2 



+ - \x\ r \/xe[-R,R]. 



In order to get a lower bound, we begin with the easier case m = 3. Recall that, when x 7^ 0, 
we use the convention x r := |x| r_1 x. First of all note that 



also, 



r x 1 "' 1 dv 1 

$r,m(x)/x r > x r =~ Vxe(0,i2]; (5.16) 



if, n/ r /o'V 1 1 jaj 1 

<&r,z{x) x = j-j- = — - > VxG -1,0 . 5.17) 

|a+ r r + 1 r(r + 1) 



So we are left to study the minimum of <& r ^(x)/x r as x varies in [-R, —1). We have: 

2/-% - 1| dy _ Ji yr~\l -y)dy + jf 1 y<-\y - 1) dy _ 



/ r (x) = $ ri3 (s)/a: r 



m^_^+ 

^ r+1 r 2 X 1 



= r(r+i) i,: 1 — -—1^+^-- vxe[-^,-i). 

\x\ r r(r +1) r + 1 r 

Notice that 

d ( 2 _ r s 1\ 2 _ r _ x 1 



s"' + - - = -s~'- x + 



ds \r(r + 1) r+1 r) r+1 r+1 

1 

whose zero is attained at sq = 2 r + 1 , so that 



f r {x) > - (2^ - 1) Vx G [— i?, —1) • (5.18) 



We finally observe that 



1 / -L- \ log 2 , 

- 2- + 1 - 1 > y 7T" • 5.19 

r V / ~ r(r + 1) v y 

By collecting together (|5.16p . (|5.17p . (|5 . 18)) and (|5.19p one gets the lower bound in (|5.15p for 
m = 3. 

Now, let us consider the case m > 3. Since the function g(s) = |s|^ m_1 - > / 2 is convex, from 
Jensen's inequality we have: 

771 — 1 



So y r x dy 



and we can bound the right hand side from below just by applying to it the above estimates 
for m = 3, so as to obtain 

%^>^ V* €[-*,()) 

|x| r ~ r 1 + 22i 2 li 

for some constant Ci = Ci(m), which together with (j5. 16|) gives the claimed lower bound in 
(|5.15p . We remain to deal with the case m G (1,3), that is when the function g(s) is no more 
convex. Again, since (|5.16p holds for all m > 1, we can restrict ourselves to analyse the ratio 
l^m^)!/!^! 7 * for x G [— R, 0). Since (m — l)/2 = a G (0, 1), straightforward calculations show 
that 

| y _l|«> |y Q -l| Vj,e[0,J?]; 
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therefore, 

\^ r (x)\ = S^y^y-i\^d y> ^yr-^-i\dy yxe[ _ m (5 20) 

\x\ r \x\ r \x\ r 

From now on, thanks to (|5.20p . by means of computations analogous to the ones developed in 
the case m = 3 one can prove that also the right hand side of (|5.20p is bounded from below by 
a constant times 1/r 2 . □ 

The next lemma is a crucial one. Basically it is a refinement of the already mentioned Lemma 
5.8 of [13] in the particular case of the functions $ r ,m(')) an d so its proof is similar to the one 
of the latter. 

Lemma 5.9. Suppose that there exists a constant Cps such that the Poincare-Sobolev inequality 

II«-^II2^<Cp5||V^|| 2 .^ WgW 1 ' 2 ^;*/,//) (5.21) 

holds for some a > 1. Then for any m > 1 and R > 1 there exists a constant = 
C*(m, v, fj,, R) such that 

||*r,m(OII 2 „.„ < || V^OH^ (5.22) 

Vr > 1/2 , V£ : U\\oo <R,S = 0, $r,m(0 G W 1 ' 2 ^; is, fi) . 

Proof. As in |13|. Lem. 5.8] we proceed by contradiction. That is, if such a constant C* did 
not exist then there would be a sequence of numbers r n > 1/2 and a corresponding sequence 
of non-identically-zero functions such that r n — > oo, Hindoo < R, £n = 0, $r„,m(£rc) G 

W 1,2 (Q',i>, fi) and 

||V* rBim (^)|| 2;/1 < \ ||*r„,m(£«)|| 2(T .„ • (5.23) 

The fact that the sequence {r n } cannot accumulate at some fixed r and therefore has to go to 
infinity is a straightforward consequence of the method of proof of |13^ Lem. 5.8]. Now, let us 
set 

a n = 1 1 $r„,m(£n) || 2a;u 



and 



From ()5.23p we get that 



l*nlL ;i , = l> l|V*n|| 2 ; M <- 



Applying the Poincare-Sobolev inequality (|5.2ip to ^> n we obtain 

pn-^n\\ 2 ^<C PS ^ (5.24) 

in particular, since H^nl^o-i/ = 1> the sequence of real numbers {^n} is bounded and therefore 
convergent to some constant cq (up to a subsequence which we do not relabel). Hence by (I5.24p 
also {^ n } converges in L 2f7 (f2; v) to such constant cq: in particular cq ^ since H^nlho-;^ = 1- 
Consider the sequence of functions 

in 



First of all, we want to prove that {2 n } converges at least pointwise to a non-zero constant. 
Indeed, again up to subsequences, we know that {^/ n } converges pointwise to the non-zero 
constant cq. Let us rewrite Z n in the following way: 



l^r„,m(Cn(x))| 



i 
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Thanks to estimate ()5. 15[) and to the fact that |£ n (x)| < R, we have: 



Co 



i 

r n 



|*n(x)|^ < l^n(x)| < 



l+lv[2S=i] 



l*n(x) 



Therefore by letting n — > oo we deduce that {|iJ n |} converges pointwise to 1, and so {Z n } 
converges pointwise to 1 if cq > or to — 1 if cq < 0, in any case to a non-zero constant. Finally 
we prove that {Z n } also converges in to such non-zero constant. In order to do that, 

thanks to Egoroff's Theorem it is enough to show that f E \Z n \ dv goes to zero as n — > oo and 
\E\ -> 0. In fact, 



I \Z n \du = [ 
Je Je 



16. 



^dz. < 



E r n 



1+1V 



Ci 



1 

'■)] 



|*n| r ™ dv < \E\ 



1 

rn 



1+1V 



r ?1 



Ci 



i 

IvD IT" 



and the assertion follows since {^ n } converges in L 2cr (Q;v), so in particular is bounded in 
L 1 (0; v). Thus {Z n } is a sequence of zero-mean functions which converge in L 1 (0; v) to a non- 
zero constant, and this is clearly absurd. Therefore the constant C* in (|5.22[) must exist. □ 



A crucial point in Lemma [5.91 lies in the fact that the constant C* does not depend on r, and 
this is necessary in order to prove the following theorem, which improves considerably the main 
result of |18| , at least in all the cases dealt with therein for which a weighted Sobolev inequality 
holds. Indeed, in |18| local uniform convergence is shown for the specific weights considered 
in the one-dimensional setting. For such weights a Sobolev-type inequality may or may not 
hold, and hence one does not expect uniform convergence in general. Our contribution is the 
identification of the role of a weighted Sobolev inequality in this connection and the functional 
analytic setting in which we prove our results. 



Theorem 5.10. Let < oo and let the Sobolev-type inequality (|1.3p hold true for some 

a > 1. Then any solution u of (|1.2p corresponding to an initial datum uq € L <?0 (O;zv), with 
qo € [1, oo) and uq = u ^ 0, converges uniformly to its mean value u. 

Proof. We shall proceed through a Moser iterative technique. Indeed, for a given g > 1, after 



multiplying equation (|1.2p by (u 



_1 and (formally) integrating by parts, straightforward 



computations lead to the following differential equation for the quantity ||w(s)||f ;i ,, where we 
set w to be the relative error w(-, s) = u(-, s)/u — 1: 



d^ 



wis 



\Q 



-q(q 



l)m\u\ m - 1 



V$#, m (w(s)) 



2:/( 



(5.25) 



Since u(r) € L°°{Q) for any given r > (recall the regularizing effect from Theorem I5.2p and 
||^(s)||oo is non-increasing (consequence of (|5.25p itself), of course there exists a suitable R > 1 
such that \w(-, s)| < R for all s > r; notice that the constant C* in Lemma [5 .91 will then depend 
on the solution u only through [|w(r)[|oo. Also, by definition, w(s) = 0. Therefore we can apply 
to the right hand side of (|5,25p the functional inequality (|5.22H (with r = g/2), which entails 

d „ , s ,,„ . p(p — l)m|u| m 
ds 



\<*)\\% v < 



2 ' 



(w(s)) 



law 



(5.26) 



Thanks to the estimate (j5. 15|) (again with r = g/2), from (|5.26p we get 
d „ . g(g - l)m\u\ m - 1 C\ / yi 2 Q 



ds 



Hs)\\^< 



e x2+2V(m-l) 



(f) 



w(s) 



< 



2a:u 



,2V(m-l) 



\w(s)f 



ag;v ' 



(5.27) 

where Q > is a suitable constant depending on m, C\, C*, \u\, but independent of g. Now, 
for a given t > r, let us set t n = r + (1 — 2~ n )(t — r) and p n +i = o~p n . Clearly, p n = cr n po 
(let pq be any real number belonging to [1, oo)), so that we have t n — > t and p n — >• oo as 



POROUS MEDIA EQUATIONS WITH HOMOGENEOUS NEUMANN BOUNDARY CONDITIONS 



21 



n — > oo. Integrating (|5.27p between t n and i n +i with g = p n and exploiting the non-expansivity 
of || Pn+1 ;jy we obtain: 

Pn 

which reads 

/ 2V(m-l) „n+l \ ST j 

lk(Wi)ll Pn+i; , < ( g 1 ("ST |^(t re )|| pre;i/ . (5.28) 

Since p n = o~ n po, for a suitable constant .K" = K(pq, m, a, Q) we can rewrite (|5.28|) as 



n+l 



that is 



ll^(*n+l)llp n+i;l /<^ p " t Pn IK*n)ll p „ ;l , , 

k(Wi)ll Pn+i; , < K^=°^ r E *= ^ |Kr)|| po;i , . 



Hence, 

IK^Hoo =Hm Ht)^;, < Ifeigf < ^T^) \\w(r)\\ mv , (5.29) 

being iiT' another suitable constant depending K, po, a and therefore on po, m, $7, v, (i, u, 
\\ w { T )\\oo- The assertion then follows by letting t — > oo in ()5.29p . □ 

The above result can be refined by showing the following analogue of Theorem 14.31 



Theorem 5.11. Let i/(fi) < oo and let the Sobolev-type inequality (|1.3p hold true for some 
a > 1. For any solution u of (|1.2p corresponding to an initial datum uq G L 1 (f2;zv) with 
uq = u 7^ i/iere exists a constant G = G(m,$l,u, fj,,u,\\uo\\i-u) > such that the following 
estimate holds: 

||«(t) -«|loo <Ge ^ Vi > 1, 

being Cp the smallest constant such that (jl.5p holds. 

Proof. One can proceed along the lines of the proof of Theorem 14.31 In order to get the expected 
exponential decay of the quantity \\u(t) —u\\2;v, one first uses estimate (|5.29p instead of (|4,12p . 
Indeed, note that the only relevant point of (|4.12p is that its right hand side goes to zero and it 
is integrable as t — > oo: both of these properties can be achieved by the right hand side of (|5.29p 
upon choosing po sufficiently next to 1 (or even pq = 1). Finally, the decay rate of the L°° 
norm is the same as for the L? norm thanks again to (|5.29p evaluated (for example) between t 
and t = t — 1/2 with po = 2. The dependence of the multiplicative constant G on the stated 
quantities follows from the constant K' in (|5.29p and from the regularity estimate (|5.2p . □ 

We remark that, obviously, Theorem 15.111 also applies to the non-weighted case. In fact the 
conclusion is the same as the one of Theorem 14.31 in terms of time decay rate. Nonetheless, in 
Theorem 14.31 the multiplicative constant G depends in a nontrivial way on the solution itself, 
while the proof of Theorem 15.111 shows that, basically, it depends on the solution only in terms 
of the initial datum (though through the constant C* of Lemma I5.9[ which is unknown) . 

We conclude giving an immediate corollary of Theorem 15.101 concerning the evolution of the 
support of (non-zero mean) solutions to (|1.2p . This topic has been widely investigated in the 
literature, see e.g. [El [25j [H [19j [28] . 

Corollary 5.12. Let u(fl) < oo and let the Sobolev-type inequality (jl.3p hold true for some 
a > 1. Then the support of any solution u(-,t) of (jl.2p corresponding to a compactly supported 
initial datum uq having non-zero mean becomes the whole f2 for all t great enough. 



We recall that in [13] . where the weighted equation (|1.2p was analysed assuming only the va- 
lidity of the Poincare inequality (|1.5p , the authors provided weights for which the corresponding 
solutions to (|1.2p do not converge uniformly to their mean value (but they do in L e {£l\v) for 
all q G [1, oo) as a consequence of (|l,5p itself). This should indicate that uniform convergence 
is strongly linked with the validity of (|1.3p . 
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5.1. Examples of weighted Sobolev inequalities. In the following we list some basic ex- 
amples of domains £1 and couples of weights (p u , p^) for which the Sobolev-type inequality (jl.3j) 
holds for appropriate values of the parameter a > 1, given below. 

• Intervals: 

o (x a 1 xf ) ) on (0,6) (b > 0): p < 1, a > -1 and all a > 1 OR P > 1, a > p - 2 and 
*e(l,f±i]; 

o xP) on (a,+oo) (a > 0): P > 1, a < -1 and all a > 1 OR p < 1, a < P - 2 
and <7 g (l, «±i ; 

o (i|logx|°,a;|logx| /3 ) on (0, c) (c G (0,1)): P > 1, a < -1 and all cr > 1 OR /3 < 1, 
q < /3 - 2 and cr G fl, |±± 

o ( e aN ;e /3M) on R: /9 > 0, a < and all a > 1 OR P < 0, a < P and cr e (l, | . 

• Bounded Lipschitz domains (N > 1): 

OR P > 1, a > P - 2 and 



JV a+N 



AT-2' JV-1 

(<5 denotes the distance function from 90). 



o (5 a ,5P): P < 1, a > -1 and a G M,min 

a G ( 1 min f a+JV " 

a t yi, mm ^ JV _ 2 , /3+N _ 2 

The Euclidean space W N (N > 1): 

o ((1 + |x|) Q , (1 + |x|)^): p > 2 - TV, a < -N and a G (l * 



N-2 



OR P < 2 — N, 



N 



a+N 



a < P — 2 and cr G ( 1 , min ^ , /3+jV _ 2 
o (e a l x l,e^l x l): p > 0, a < and a G (l, 
1, min 



AT-2 



OR /3 < 0, a < P and a G 



AT a 



7V-2' /3 y 

The above examples can be obtained by applying Theorem 1.4 of [7] (which gives necessary 
and sufficient conditions on the weights so that (j 1 . 3 1) holds) in the one-dimensional case and 
by exploiting the results of |21| in the ^-dimensional case, see in particular Chapter 19 for 
bounded Lipschitz domains and Chapters 20, 21 for unbounded domains (the validity of (j 1 . 3 f) 
here is a consequence of the compact embeddings discussed therein). 
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